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Abstract — We study the performance of binary spatially- 
coupled low-density parity-check codes (SC-LDPC) when used 
with bit-interleaved coded-modulation (BICM) schemes. This 
paper considers the cases when transmission takes place over 
additive white Gaussian noise (AWGN) channels and Rayleigh 
fast-fading channels. The technique of upper bounding the 
maximum-a-posteriori (MAP) decoding performance of LDPC 
codes using an area theorem is extended for BICM schemes. 
The upper bound is computed for both the optimal MAP 
demapper and the suboptimal max-log-MAP (MLM) demapper. 
It is observed that this bound approaches the noise threshold of 
BICM channels for regular LDPC codes with large degrees. The 
rest of the paper extends these techniques to SC-LDPC codes 
and the phenomenon of threshold saturation is demonstrated 
numerically. Based on numerical evidence, we conjecture that the 
belief-propagation (BP) decoding threshold of SC-LDPC codes 
approaches the MAP decoding threshold of the underlying LDPC 
ensemble on BICM channels. Numerical results also show that 
SC-LDPC codes approach the BICM capacity over different 
channels and modulation schemes. 

Index Terms — BICM, Rayleigh fast-fading, density evolution, 
GEXIT curves, LDPC codes. 



I. Introduction 

The phenomenon of threshold saturation was introduced 
by Kudekar et al. [ 1 ] to explain the impressive performance 
of convolutional low-density parity-check (LDPC) ensembles 
[2), 0. These codes are essentially terminated convolutional 
codes with large memory, which admit a sparse parity-check 
matrix representation. One way to construct these codes is 
to "spatially-couple" an underlying LDPC ensemble, resulting 
in a spatially-coupled LDPC (SC-LDPC) ensemble. It was ob- 
served that the belief -propagation (BP) threshold of a spatially- 
coupled ensemble is very close to the maximum-a-posteriori 
(MAP) threshold of its underlying ensemble; a similar state- 
ment was formulated independently, as a conjecture in J4). 
This phenomenon has since been called "threshold saturation 
via spatial coupling". Kudekar et al. prove in Qj that threshold 
saturation occurs for the binary erasure channel (BEC) and a 
particular class of underlying regular LDPC ensembles. For 
general binary-input memoryless symmetric (BMS) channels, 
threshold saturation was empirically observed first 0, [Q and 
then analytically shown Q, 0. It is known that the MAP 
threshold of regular LDPC codes approaches the Shannon 
limit for binary memoryless symmetric (BMS) channels with 
increasing left degree, while keeping the rate fixed (though 
such codes have a vanishing BP threshold) (TJ. So, spatial 
coupling provides us with a technique to construct a single 
capacity approaching code ensemble for all BMS channels 



with a given capacity. This technique is indeed very general 
and has since been applied to a broad class of graphical 
models. A good summary of recent applications of spatial 
coupling can be found in |8|. 

In this paper, we evaluate the performance of spatially- 
coupled LDPC codes using BICM schemes for transmission 
over additive white Gaussian noise (AWGN) and Rayleigh 
fast-fading channels. The noise threshold, a.k.a. the Shannon 
limit, for bit-interleaved coded-modulation (BICM) schemes 
can be computed using Monte-Carlo simulations via the gener- 
alized mutual information (GMI) [9|. This method can be used 
to compute the information theoretic limits for different subop- 
timal BICM schemes and is briefly reviewed in Section Hl-BI 
We review density evolution (DE) for BICM schemes, de- 
scribed in JTU] Sec. 5.2], in Section III-CI We note that the 
above DE can be greatly simplified by using the Gaussian 
mixture approximation for the BICM bit-channels presented 
in ifTTIl . to obtain approximate thresholds. Section UlTl extends 
the GEXIT analysis and the upper bounding technique on 
the MAP decoding threshold for BICM schemes. Section HVl 
extends the analysis to SC-LDPC codes. The DE results of SC- 
LDPC codes are presented in Section [V] and some concluding 
remarks are given in Section [VI] 

II. Background 

A. The BICM Model 

BICM is a practical approach to coded modulation and was 
introduced by Zehavi in lfl2l . A comprehensive analysis for 
BICM is provided in iflOll . which is an excellent reference 
for BICM. We now briefly describe the BICM model and 
the problem setup. Consider transmission over a memory-less 
channel with input alphabet X (with \X\ =2 M ,M £ N) and 
output alphabet C. We use uppercase letters (e.g., X, Y) to 
denote random variables and lowercase letters (e.g. x, y) to 
denote their corresponding realizations. The channel output is 
given by 



Y = AX + Z, 



(1) 



where X 6 X, Y G C, and Z is additive Gaussian noise with 
variance a 2 i.e., Z ~ CAf(0, er 2 ). We consider the cases of 
no fading (A = 1) and Rayleigh fast-fading (A - CAA(0, 1)). 
Furthermore, we assume that the receiver has perfect channel 
state information for simplicity. The analysis can be easily 
extended to the case when the receiver does not have access 
to the channel state information IT3l . Ifl4l Sec. 5.1]. 
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Fig. 1, The BICM system model. The codeword is mapped to a symbol 
X £ X using the mapper. The channel output Y is then passed through a 
demapper which performs the symbol-to-bit metric calculation. The decoder 
output can be optionally fed back to the demapper as apriori information for 
BICM-ID. 



A Bernoulli-(l/2) source is encoded using an LDPC code 
chosen uniformly at random from the standard irregular en- 
semble LDPC(N,A,p) [14, Ch. 3]. Here, \{x) = 'Zi^' 1 
is the degree distribution (from an edge perspective) cor- 
responding to the variable nodes and p(x) = Zi Pi- 1-1 
is the degree distribution (from an edge perspective) of the 
parity-check nodes in the decoding graphjj The coefficient Xi 
(resp. pi) gives the fraction of edges that connect to variable 
nodes (resp. parity-check nodes) of degree i. Likewise, let Li 
be the fraction of variable nodes with degree i and define 
L(x) — J2i LiX 1 . The design rate of the LDPC code is given 
by 



i(x, P )=i- f j p{x)dx . 

Jo X(x)dx 



The blocklength TV is assumed to be a multiple of M, where 
groups of M bits are mapped to a symbol in X and then 
transmitted over the channel. At the receiver, a demapper 
first performs the symbol-to-bit metric calculation based on 
the received symbol Y, and the metrics are then passed to 
the decoder. One can also perform the symbol-to-bit metric 
calculation iteratively, by using the decoder output as apriori 
information at the demapper. This scheme is commonly known 
in the literature as BICM iterative detection (BICM-ID). The 
block diagram of a general BICM system is shown in Fig. Q] 
A BICM scheme is specified by the bit-to-symbol mapper 
and the demapper. Throughout this work, we consider square 
quadrature amplitude modulation (QAM) constellations with 
the Gray mapping scheme, and the optimal MAP and subop- 
timal max-log-MAP (MLM) demappers. 

B. Noise Threshold of BICM Channels 

Consider the case when the demapper calculation is not 
updated between iterations. The performance of the BICM 
scheme (for optimal demappers) is given by the capacity of a 
set of parallel independent channels |fT31 . It was also character- 
ized in terms of the generalized mutual information (GMI) by 
viewing the BICM decoder as a mismatched decoder [ 16 1. The 
GMI analysis was used in ifTTIl to compute the performance 
of BICM in the presence of suboptimal demappers. The 
achievable information rate of a given BICM scheme can be 
computed using Monte-Carlo simulations via the GMI |9). 
Consider a BICM channel with M -bits per symbol. The I- 
curve was introduced in [9] and can be computed for the m-th 

'The edges of the variable nodes connected to the demapper are not 
included in the degree profile. 
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Fig. 2. The Tanner graph at the decoder for BICM channels. The red 
diamonds represent the demapper nodes, the blue circles and black squares 
represent the variable and check nodes respectively. Each demapper node is 
connected to M variable nodes. 



bit level via 

I m (s; a) = 1- E x , r log(l+exp((26 m (X) - l)K m (Y(a))s) ) , 

where X and Y are the channel input and output respectively 
(we have used the notation Y(a) to make explicit the depen- 
dence on the channel noise variance), b m (X) is the m-th bit 
label of symbol X and A m (Y) is the log-likelihood ratio of the 
m-th bit of the symbol after passing through the demapper |0. 
The /-curve of the BICM channel is then computed as 

M-l 

I(s:a) = 2J I m (s;a). 

The achievable information rate of the BICM scheme is equal 
to the GMI, given by 1(a) — max s I(s; a). This enables us 
to compute the achievable information rate of BICM schemes, 
for different modulation schemes, bit-to-symbol mappings and 
demappers. The noise threshold for error-free transmission at a 
given transmission rate R can then be computed for a specific 
BICM scheme by a* = I~ X (R). 

For the case of BICM-ID, the capacity would be equal to 
the coded modulation capacity, when the input alphabet is 
restricted to X. The noise threshold can be computed for this 
case similarly. Let 1(a) — I(X;Y(a)), when the input is 
uniformly distributed. Then the noise threshold for BICM-ID 
is given by a* — I~ 1 (R). 



C. Density Evolution 

We begin this section by first introducing some notation. 
Let Vi, Cj and dk denote the variable, check and demapper 
nodes respectively. Let n(k, m) = (k — 1)M + m = i, be the 
mapping from the demapper nodes to the variable nodes i.e., 
the m-th bit of demapper node k is connected to variable node 
i. When the symbol index k is understood from context, we 
write 7r(m) = i, and m = 7r~ 1 (i). The m-th bit corresponding 
to x G X is denoted by b m (x), and the set of symbols where 
the m-th bit is zero (one) is denoted by X™ (X" 1 ). 

The factor graph structure at the joint decoder is shown in 
Fig. |2 The joint decoder proceeds by performing one round of 
decoding for the LDPC code followed by a demapper update. 
This is the schedule for BICM-ID. To reduce the complexity 
non-iterative detection is used, where the demapper update is 
not performed between iterations. 

Let pif^, Hcj^n, Pv^d k and Pdl^ be the messages 
from the bit node to check node, check node to variable node, 
variable node to demapper node and demapper node to variable 
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node during iteration I respectively. All the messages are in 
log-likelihood ratio domain. The message passing rules at the 
variable and check nodes are the standard rules and their 
description is omitted. Using the notation di to denote the 
set of check nodes connected to variable node i, the message 
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Let to be the bit index corresponding to variable node Vi i.e., 
to = 7r _1 (i) and B m be the random variable corresponding 
to that bit. The bit probabilities in the £-th iteration can be 
computed using the variable node to demapper node messages 
via 



Pr(B m = 0) 



A - 



-, Pr(K 



-I) 



1 



So, the demapper to variable node message is given by 



CO 



log 



£* e *r p(vk\ x ) Uiyim Pr ( B 7r(o =K(i) {x)) ' 

(2) 



where B n m is used to denote the bit corresponding to variable 
node Uwfe,;) and is an abuse of notation. The above message 
passing rule is for the optimal MAP demapper. The rule for 
the MLM demapper is obtained by performing the standard 
approximation of the above equation to reduce complexity. The 
variable nodes % = 1, ■ ■ ■ , N can be grouped into equivalence 
classes via the function n(k,m) i.e., let V m = {i>i|7r(fc, to) = 
i, k = 1, • • • ,N/M} denote the set of all variable nodes 
connected to the TO-th bit of the demapper nodes. Denote the 
density of messages emanating from the variable nodes in V m 
to the check nodes at iteration I by a™, conditioned on the 
transmission of an all-zero codeword. Note that the all-zero 
codeword assumption is not valid, but we can still use DE 
with standard symmetrizing techniques 03], |[T8] Ch. 7]. The 
transformation of the densities of the incoming messages at 
the check node and variable node are denoted by EH and © 
respectively (see discussion in Ifl4l p. 181]). For a density x, 
we denote 

x mn = X ffl X ffl • • • ffl X, 



and likewise for x®". Using this notation, define A(x) = 



i>(i-l) 



p(x) = E^x ffl(i - 1} and L(x) = 



The density of messages at the input to the check nodes is 
given by 



M 

— T 



(3) 



We call the average density of messages from the variable 
node to check node at iteration I. The density of messages 
from the variable node to the demapper node is then given by 
L(p(a^')). Let m (-;<x) be the demapper density transforma- 
tion operator of the rn-th bit corresponding to (O. Then, the 



density evolution equations are given by 

a r ^ 1 )=0 m (L(p(aW));a)©A( /5 (aW)) 



1 



M 



M ^ m ' 

m— 1 



(4) 



from which one obtains the recursion 

a(^^(LWaW));,)iA(p(aW)), 

where </>(■; a) maps the incoming density at the demapper node 
to the average output density i.e., 

1 M 

m— 1 

One can use the 'M' equations to perform DE for 
protograph based LDPC codes to design bit mappings for 
optimal performance, similar to [19 | where the authors use 
PEXIT curves for optimization. This function does not have 
a closed form expression and can be computed using Monte- 
Carlo simulations. 

III. GEXIT Curves for BICM 

In this section, we derive an expression for the BP-GEXIT 
curve for LDPC codes for BICM schemes. Using the BP- 
GEXIT curve and the area theorem, an upper bound is derived 
for the MAP decoding threshold of LDPC codes for BICM 
schemes. As defined in the previous section, let 7r(fc, m) = 
(k — 1)M + to = i, be the mapping from the demapper 
nodes to the variable nodes. In this section, boldface uppercase 
letters (e.g. X,Y) are used to denote random vectors and 
X^k to denote the vector with all elements of X except the 
fc-th element. Let x[ m ' be the m-th bit of symbol k and let 
i = 7r(fc, to). Throughout this section, variable node V{ shall be 
denoted by y i a tne function n. Consider transmission over 
the BICM channel family (Q]) parametrized by the normalized 
channel entropy per bit, given by 



a 



We note that a G [0,1] and that the channel family is 
complete and degraded with respect to a. Let X 6 %n/m 
be the transmitted vector and Y be the output of the channel. 
Following the definition for BMS channels |[T4l Ch. 4], the 
GEXIT function for BICM channels is defined as 



1 AH(X\Y{a)) 



N da 
and satisfies an area theorem by definition: 



1 



N 
R(A,p) 



(H(X\Y(l))-H(X\Y(0))) 



(5) 



It is convenient to assume that symbol k is transmitted 
through a channel with parameter a^, and that each is 
further characterized by a common parameter a in a smooth 
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and differentiable manner. For the case under consideration, 
we simply have a k = a. Then define the k-th GEXIT function 

dH(X\Y(a u --- ,a N/M )) 



gfc(ai 



,ckn/m) = 



dai 



So, the GEXIT function is given by 

N/M 



(«) = —>. Ek{ai,--' ,a N/M ) — 



N ^ 

k=l 



da k 
da 



Lemma 1: Consider transmission using an LDPC code from 
the ensemble LDPC(A, A,p) over the BICM channel with 
parameter a. Define (j> k {y^ k ) = {Px k \Y„ k (%\y~k),x e X] 
and let & k (Y^ k ) be the corresponding random variable. Then, 
the fc-th GEXIT function is given by 



gfc(a) = G(x fc (u);a) = ^ p{x k ) / x Xk (u)n Xk (u)di 

... , T.' •' U 



x k £X 



(6) 



where x fc (w) = {x Xfe (u), x k e -Y}, and x X(c (tt) = ?>(^>fe|iKfe) is 
the distribution of <f) k assuming that X k = x k was transmitted, 
and the GEXIT kernel is given by 

T, x >£x u [ x ']p(y\ x ') . 



K x(u) = J —p(y\x)\og 2 



-dy, (7) 



t y — tt[x]p(y|x) 

where u[j] denotes the j-th component of it. 

Proof: It can be verified that is the extrinsic MAP 
estimator of X k . We have 

. dH(X\Y(a u --- ,a n )) dH(X k \Y k) $ k ) 

Sfc( a ) = a = n ■ 

da k da k 

For notational convenience, we omit the dependence of X k , 
Y k and $ k on the symbol index k whenever possible. The 
conditional entropy of X k is given by 

v(x,y,4>) 



H{X\Y^) = - /j>(x ;2 /,0)log. 



dydcj) 



:J2p(x)I p (<p\ X )\ / P (,ix)iog 2 E -'^gg^ |a:/) d,j 



xEX 



p(x\<p)p(y\x) 



This follows by noting that p(x k ,y k ,(f> k ) = 
p{y k \x k )p((f> k \x k )p(x k ). The result now follows by noting 
that p(x fc |0 fc ) =_p(x fc |y^fc). 

Assuming that each bit in the symbol x k is independent 
(which is true asymptotically as iV -> oo), we have 

f M ~\ 



\y~ k ),x k e X) 



{fx k (v!, ■ ■ ■ ,V M ),X k e x}, 



where v n 



log 



p(x [ k m] =+l\y^ k ) 



From dHJ, we see that u is 



' p( x k , =-i|»~k) 

completely characterized by Ui, • • ■ , dm. So, Q7p is henceforth 
interpreted in terms of the log-likelihood ratios v m . ■ 
The density x k (u) in (|6]l is hard to compute, so instead 
one can use the BP estimate to compute the density for the 
asymptotic limit N — > oo. The curve obtained by using the BP 
estimate is called the BP-GEXIT function g BP (a). Let x (£) (w) 



denote the density of the log-likelihood ratio, conditioned on 
the transmission of the all-zero codeword, emitted from the 
variable nodes to the detector nodes, during iteration i. For 
a fixed £, as the blocklength — s- oo, we have p x (v m ) = 
If F[x] is the density transformation operator 
corresponding to the map p n- e M /(l + e M ), then we can write 



x(ix) = {p(ix|x),x e X} 



= jll F[x](x^v m ),xeX^ 



(9) 



where u is given by (0. The BP-GEXIT function g BP 
is computed as follows: For a given channel parameter a, 
compute the fixed point of density evolution, say a. Then, 

g BP ( a ) = G(F[L(p(a))];a). 

One can now calculate an upper bound on the MAP decod- 
ing threshold. The following procedure is now fairly standard 
and the details can be found in lfl4l Sec. 4.12]. It can be shown 
that the GEXIT functional G preserves degradation. So, by the 
optimality of the MAP decoder, the GEXIT function always 
lies below the BP-GEXIT function i.e., g(a) < g BP (a). Let 
a be the largest positive number such that 



,bp 



(a)da — R(A, p). 



From the properties of the BP-GEXIT function and the GEXIT 
function, we have a MAP < a. Following Q, we refer to a 
as the area threshold. 

In this work we compute the BP-GEXIT curves for the case 
when there is no demapper update between iterations. The BP- 
GEXIT functions and the area threshold, for different BICM 
schemes are shown in Figures H] |5] |6] The parameter a was 
chosen to be the normalized channel entropy per bit. 

IV. Spatially-Coupled LDPC Codes 

In this section, we describe the spatially-coupled 
(l, r,L, w) ensemble introduced in (TJ. The variable 
nodes are placed at positions [— L, L] and the check nodes 
are placed at positions [— L, L + w — 1] (w can be thought of 
as a "smoothing" parameter). Each of the 1 connections, of a 
variable node at position i, are uniformly and independently 
chosen from [i, i + w — 1] as shown in Fig. [3] There are two 
consequences of this coupling - threshold increase and rate 
loss [1]. The design rate of a spatially-coupled (l,r,L, w) 
ensemble is given by 

l W + l-2Eto(^) r 
r 2L + 1 



R(l,r, L,w) = 1- 



(8) Let be the average density, in the spirit of (|3), emitted by 



the variable nodes at position i. Set a^ = A +QO , the Dirac 
delta function at infinity, for i £ [— L,L]. The DE for this 
ensemble can be written as 



^ W — 1 W — 1 



(10) 



3=0 



k =0 
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i — w + 1 i — 1 i i + 1 i + w — 1 




— 1 iu — 1 



Fig. 3. A portion of a generic spatially-coupled code. The variable nodes 
are shown as blue rilled circles and the check nodes are respresented by 
black squares. The variable node at position i is connected to check nodes at 
positions [i, i+w — 1]. The permutations 7r; are chosen uniformly at random. 
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Fig. 4. The BP-GEXIT curve of a (3, 6) code for QPSK modulation and 
the upper bound on the MAP threshold computed via the area theorem. Also 
shown is the BP-GEXIT curve of the (3, 6, 32, 4) spatially-coupled ensemble. 
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Here L(x) = x®\A(x) = x® 1 " 1 , and p(x) = 
a channel parameter a, let a = [a_£, • • • , a^] denote the 
fixed point implied by ( fTOb . Using the technique developed in 
Section [HI] one can compute the GEXIT curves for spatially- 
coupled LDPC codes from © along the lines of Define 
the GEXIT functional for spatially-coupled codes via 



1 L 

G(a) = V G(ai), 

' 2L + 1 ^ K h 

i=-L 



where G(a) is defined in ©. The BP-GEXIT function 
of spatially-coupled codes is given in parametric form by 
(a,G(I)). 

V. Results 

We consider the case when there is no demapper update 
between iterations. The BP-GEXIT curve of a (3,6,64,4) 
spatially-coupled ensemble is shown in Figures [4] [5] |6] for 
different BICM schemes. Figures [4] [5] |6] also demonstrate the 
threshold saturation phenomenon, wherein the BP threshold of 
the spatially-coupled ensemble is close to an intrinsic threshold 
of the underlying ensemble. For the optimal MAP demapper, 
this intrinsic threshold is the area threshold computed using 
the BP-GEXIT curve of the underlying ensemble. From Fig. [6] 
we observe that for the suboptimal MLM demapper, the 
BP threshold of the spatially-coupled code crosses the area 
threshold of the underlying ensemble. 

It is known that increasing the left degree of regular LDPC 
codes, while keeping the rate constant, pushes the MAP 
threshold of the ensemble towards the noise threshold for BMS 
channels. Based on this, we compute the BP thresholds of a 
(4,8,64,4) and (6,12,64,4) spatially-coupled ensemble and 
as seen from Tables U and HH the gap to the BICM noise 
threshold indeed becomes smaller. The BP thresholds and 
the noise thresholds for the Rayleigh fast fading channel are 
shown in Tables HI and Hll The tables also show the asymptotic 
gap as the rate loss tends to zero, i.e., when L — s> oo. 
Based on this numerical evidence, we conjecture that the 
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Fig. 5. The BP-GEXIT curve of a (3, 6) code for 16QAM modulation with 
a MAP detector and the upper bound on the MAP threshold computed via 
the area theorem. Also shown is the BP-GEXIT curve of the (3, 6, 32, 4) 
spatially-coupled ensemble. 



ensemble of spatially-coupled LDPC codes with large left 
degrees universally approach the noise threshold for different 
BICM schemes. 

The decoding thresholds of coding schemes for BICM 
can be improved by using BICM-ID. For comparison, the 
thresholds of the (4,8,64,4) spatially-coupled ensemble is 
computed for BICM-ID. The thresholds are computed for 
64QAM modulation using a MAP demapper and we con- 
sider transmission over both AWGN and Rayleigh fast-fading 
channels. In this case, the demapper update was performed 
once for every 100 iterations of the spatially-coupled code. As 
expected, the gap to the noise threshold reduces to 0.09 dB 
from 0.20 dB and 0.24 dB for AWGN and Rayleigh fast-fading 
channels respectively. This improved performance comes at 
the expense of increased complexity. 
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TABLE I 

Performance of various spatially-coupled ensembles over BICM channels with a MAP demapper. 





BP Thresh./ Gap /Asympt. Gap {E b /N a in dB) 


Mod. / Chan. 


Noise Thresh. (E b /N in dB) 


(3,6,64, 4) 


(4, 8, 64, 4) 


(6,12,64,4) 


QPSK / AWGN 


0.17 


0.57/0.40/0.31 


0.54/0.37/0.11 


0.33/0.16/0.06 


16QAM / AWGN 


2.27 


2.71/0.44/0.35 


2.49/0.22/0.13 


2.43/0.16/0.07 


64QAM / AWGN 


4.67 


5.11/0.44/0.35 


4.87/0.20/0.11 


4.82/0.15/0.05 


QPSK / Fading 


1.83 


2.27/0.44/0.35 


2.04/0.21/0.12 


2.00/0.17/0.07 


I6QAM / Fading 


4.11 


4.56/0.45/0.36 


4.34/0.23/0.14 


4.29/0.18/0.08 


64QAM / Fading 


6.62 


7.11/0.49/0.40 


6.86/0.24/0.15 


6.80/0.18/0.08 



TABLE II 

Performance of various spatially-coupled ensembles over BICM channels with an MLM demapper. 





BP Thresh. /Gap /Asympt. Gap (E b /N n in dB) 


Mod. / Chan. 


Noise Thresh. (E b /N in dB) 


(3,6,64, 4) 


(4, 8, 64, 4) 


(6,12,64, 4) 


QPSK / AWGN 


0.17 


0.57/0.40/0.31 


0.54/0.37/0.11 


0.33/0.16/0.06 


16QAM / AWGN 


2.29 


2.70/0.41/0.32 


2.51/0.22/0.13 


2.47/0.18/0.08 


64QAM / AWGN 
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4.96/0.25/0.15 


QPSK / Fading 


1.83 
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2.00/0.17/0.07 


16QAM / Fading 


4.17 


4.63/0.46/0.37 


4.41/0.24/0.15 


4.36/0.19/0.09 


64QAM / Fading 


6.73 


7.26/0.53/0.44 


7.01/0.28/0.19 


6.95/0.22/0.12 
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Fig. 6. The BP-GEXIT curve of a (3, 6) code for 16QAM modulation with 
an MLM detector and the upper bound on the MAP threshold computed via 
the area theorem. Also shown is the BP-GEXIT curve of the (3, 6, 32, 4) 
spatially-coupled ensemble. 



VI. Concluding Remarks 

Spatially-coupled LDPC codes have shown promising re- 
sults for a large class of graphical models. In this work, 
we study their performance on BICM channels and validate 
the conjecture that the phenomenon of threshold saturation 
is indeed very general. We extend the GEXIT analysis of 
LDPC codes for BICM schemes. This enables one to bound 
the performance of the MAP decoder of LDPC codes. The 
area threshold, which upper-bounds the MAP decoding perfor- 
mance, of LDPC codes is computed for different demappers 
and modulations. Using these tools, we numerically demon- 
strate the phenomenon of threshold saturation for these chan- 
nels. We note that when using suboptimal demappers (like the 
MLM demapper), the threshold saturates towards an intrinsic 
threshold of the suboptimal demapper and the upper bound 
computed using GEXIT curves is no longer tight. This is con- 



sistent with previous results for spatially-coupled systems with 
suboptimal component decoders [20|-[22|. The performance 
also improves significantly when used with BICM-ID and the 
thresholds of SC-LDPC codes approach the noise threshold 
with smaller degrees. These asymptotic results demonstrate 
that SC-LDPC codes approach the noise threshold of different 
BICM schemes. 

Irregular LDPC codes and multi-edge type LDPC codes 
were designed for BICM schemes in |18|, [23 1 with excellent 
thresholds. In [19|, the authors design families of protograph- 
based LDPC codes for BICM schemes with excellent thresh- 
olds of around 0.2-0.4 dB from the BICM noise threshold. 
For finite L, spatially-coupled codes compare favorably with 
the previously optimized LDPC ensembles in-spite of the rate 
loss incurred due to finite L. 



References 



[i] 



[2] 



S. Kudekar, T. J. Richardson, and R. L. Urbanke, "Threshold saturation 
via spatial coupling: Why convolutional LDPC ensembles perform so 
well over the BEC," IEEE Trans. Inform. Theory, vol. 57, no. 2, pp. 
803-834,2011. 

J. Felstrom and K. S. Zigangirov, "Time-varying periodic convolutional 
codes with low-density parity-check matrix," IEEE Trans. Inform. The- 
ory, vol. 45, no. 6, pp. 2181-2191, 1999. 
[3] M. Lentmaier, A. Sridharan, K. S. Zigangirov, and D. J. Costello, "Ter- 
minated LDPC convolutional codes with thresholds close to capacity," 
in Proc. IEEE Int. Symp. Inform. Theory, Adelaide, Australia, 2005, pp. 
1372-1376. 

[4] M. Lentmaier and G. P. Fettweis, "On the thresholds of generalized 
LDPC convolutional codes based on protographs," in Proc. IEEE Int. 
Symp. Inform. Theory, Austin, TX, 2010, pp. 709-713. 
[5] M. Lentmaier, A. Sridharan, D. J. Costello, and K. S. Zigangirov, 
"Iterative decoding threshold analysis for LDPC convolutional codes," 
IEEE Trans. Inform. Theory, vol. 56, no. 10, pp. 5274-5289, Oct. 2010. 
S. Kudekar, C. Measson, T. Richardson, and R. Urbanke, "Threshold 
saturation on BMS channels via spatial coupling," in Proc. Int. Symp. 
on Turbo Codes & Iterative Inform. Proc, Sept. 2010, pp. 309-313. 
[7] S. Kudekar, T. Richardson, and R. Urbanke, "Spatially coupled ensem- 
bles universally achieve capacity under belief propagation," in Proc. 
IEEE Int. Symp. Inform. Theory, July 2012, pp. 453^457. 
[8] S. Kumar, A. Young, N. Macris, and H. D. Pfister, "A proof of threshold 
saturation for spatially-coupled LDPC codes on BMS channels," in Proc. 
50th Annual Allerton Conf. on Commun., Control, and Comp., Oct. 2012. 



[6] 



7 



[9] T. T. Nguyen and L. Lampe, "Bit-interleaved coded modulation with 
mismatched decoding metrics," IEEE Trans. Commun., vol. 59, no. 2, 
pp. 437^147, Feb. 2011. 

[10] A. Guillen i Fabregas, A. Martinez, and G. Caire, Bit-Interleaved Coded 
Modulation. Now Publishers Inc., 2008, vol. 5, no. 1-2. 

[11] A. Alvarado, L. Szczecinski, R. Feick, and L. Ahumada, "Distribution 
of L-values in gray-mapped M 2 -QAM: closed-form approximations and 
applications," IEEE Trans. Commun., vol. 57, no. 7, pp. 2071-2079, July 
2009. 

[12] E. Zehavi, "8-PSK trellis codes for a Rayleigh channel," IEEE Trans. 

Commun., vol. 40, no. 5, pp. 873-884, 1992. 
[13] J. Hou, P. H. Siegel, and L. B. Milstein, "Performance analysis and 

code optimization of low density parity-check codes on Rayleigh fading 

channels," IEEE J. Select. Areas Commun., vol. 19, no. 5, pp. 924-934, 

May 2001. 

[14] T. J. Richardson and R. L. Urbanke, Modern Coding Theory. Cambridge 
University Press, 2008. 

[15] G. Caire, G. Taricco, and E. Biglieri, "Bit-interleaved coded modula- 
tion," IEEE Trans. Inform. Theory, vol. 44, no. 3, pp. 927-946, May 
1998. 

[16] A. Martinez, A. Guillen i Fabregas, G. Caire, and F. Willems, "Bit- 
interleaved coded modulation revisited: A mismatched decoding per- 
spective," IEEE Trans. Inform. Theory, vol. 55, no. 6, pp. 2756-2765, 
June 2009. 



[17] J. Jalden, P. Fertl, and G. Matz, "On the generalized mutual information 
of BICM systems with approximate demodulation," in Proc. IEEE 
Inform. Theory Workshop, Jan. 2010, pp. 1-5. 

[18] J. Hou, P. H. Siegel, L. B. Milstein, and H. D. Pfister, "Capacity- 
approaching bandwidth-efficient coded modulation schemes based on 
low-density parity-check codes," IEEE Trans. Inform. Theory, vol. 49, 
no. 9, pp. 2141-2155, Sept. 2003. 

[19] T. V. Nguyen, A. Nosratinia, and D. Divsalar, "Threshold of protograph- 
based LDPC coded BICM for Rayleigh fading," in Proc. IEEE Global 
Telecom. Conf, Dec. 2011, pp. 1-5. 

[20] K. Takeuchi, T. Tanaka, and T. Kawabata, "Improvement of BP-based 
CDMA multiuser detection by spatial coupling," in Proc. IEEE Int. 
Symp. Inform. Theory, St. Petersburg, Russia, July 2011, pp. 1489-1493. 

[21] C. Schlegel and D. Trahachev, "Multiple access demodulation in the 
lifted signal graph with spatial coupling," in Proc. IEEE Int. Symp. 
Inform. Theory, St. Petersburg, Russia, July 2011, pp. 2989-2993. 

[22] Y.-Y. Jian, H. Pfister, and K. Narayanan, "Approaching capacity at high 
rates with iterative hard-decision decoding," in Proc. IEEE Int. Symp. 
Inform. Theory, July 2012, pp. 2696-2700. 

[23] H. Sankar, N. Sindhushayana, and K. Narayanan, "Design of low-density 
parity-check (LDPC) codes for high order constellations," in Proc. IEEE 
Global Telecom. Conf, vol. 5, Nov. 2004, pp. 3113-3117. 



